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5t , Abstract. A suggestion is made for distinguishing 2N and 6q short range 

correlations within the deuteron. The suggestion depends upon observing 
high momentum backward nucleons emerging from inelastic electromag- 
netic scattering from a deuteron target. A simple model is worked out to 
see the size of effects that may be expected. 
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1 Introduction 

In deep inelastic scattering upon deuterons or heavier nuclei, nucleons or other hadrons 
can emerge backward to the direction defined by the incoming photon in the target rest 
frame. The backward nucleons do not, we believe, come from the nucleon or the quark 
that was struck |]I|. Rather, the backward nucleons come from the debris that remains 
after the item that was struck is driven strongly forward. We have commented 0, ^ 
upon using neutrino production of backward protons [^, ^, ^, 0] to explore short distance 
quark configurations and here wish to use deep inelastic electromagnetic interactions to 
similar ends. 

We distinguish two extreme cases. In a 2-nucleon or 2N correlation the nucleons 
maintain their characters as neutrons and protons no matter what their separation is. 
To get a backward proton from the deuteron, one must strike the neutron. This will 
break up the bound state, and the proton will emerge with the Fermi momentum it has 
at the moment of breakup, and the Fermi momentum will be backwards in the case of 
interest. 

The other extreme case is the 6-quark or 6q cluster. Here we mean a "kneaded" 
6q, uuuddd, object with all the quarks in relative S-states. The fiavor-spin-color wave 
function is unique (for overall deuteron quantum numbers) and is not equivalent to two 
nucleons lying at the same spatial point. Emitting a backward proton begins with one 
quark being struck and driven forward. The proton must be formed out of the remaining 
5 quarks, plus possible higher Fock components, and the process of forming hadrons we 
refer to as the "fragmentation" of the 5-quark residuum. The term fragmentation follows 
common usage for the production of hadrons from any color non-singlet QCD object, 
quark and gluon jets being the most familiar. In the breakup of this 5-quark residuum, 
we persist with the nomenclature "fragmentation" although recombination may be the 



process chiefly at work. In any case, the 6q model can produce a backward proton spec- 
trum which agrees with data from neutrino reactions, for backward hemisphere proton 
momentum above about 300 MeV 0. However, so can the 2N model with enhanced 
high momentum components [^ . We need a more detailed indicator to test the two mod- 
els. For weak interactions, ratios of neutrino and antineutrino induced backward proton 
production cross sections cancel much of what is unknown in either the 2N or 6q model 
and yet give predicted results which are not the same for the two models [0. Similar 
opportunities should also exist for electromagnetic interactions. 

Here we suggest a characteristic in the spectrum of backward nucleons which is re- 
liably different for the two scenarios of the short range configuration, and which allows 
electromagnetic experiments to adjudicate between them. A suggestion for a ratio to ex- 
amine at fixed backward proton momentum and varying Bjorken x is given in section 2, 
along with some numerical estimates of the differences between the two models. Further 
comments on possibilities at CEBAF, Fermilab, or elsewhere are made in section 3. Inci- 
dentally, an observation regarding changes in average Bjorken x with varying backward 
nucleon momentum, which was originally made in the context of the 2N model, should 
be reasonably true for most any model and we also comment on this in section 3. We 
conclude in section 4. 

2 Signatures of short distance correlations 

2.1 General 

The cross section for inelastic electromagnetic scattering of a lepton from a stationary 
target is (in the scaling region and neglecting aL/o-r) 

da 4:71 at rriN E 



dx dy Q^ 



(l + (l-yf) F,{x,Q') (1) 



where E and E' are the incoming and outgoing lepton energies, z/ is the difference between 
them, y = u/E, x = Q"^ /2mNU, and F2 (whose Q^ dependence will generally be tacit) is 

F2{x) = Ee'a;g,(x) (2) 

where Cj is the quark charge in units of proton charge and qi is the distribution function 
for a quark of flavor i. 

2.2 Backward nucleons from 2N correlations 

We will speak of observing a backwards proton for the sake of definiteness; observing a 
backwards neutron is essentially similar. 

Some things change in the above formula when a backwards proton is observed. The 
neutron, which is the struck particle, is not at rest in the lab frame. Then the momentum 
fraction of the struck quark relative to the neutron is not x but rather ^, 

where a = {Ep+p^)/mN is the light cone momentum fraction of the backward proton with 
p^ positive for a backward proton. One should also replace itlnE by its corresponding 
Lorentz invariant, 

m]\rE -^ Pn ■ k = mj^E{2 — a). (4) 

Also, if we want a backward proton, there is a further factor of the probability density 
for finding the proton with its observed momentum at the moment the neutron was struck. 
Thus 

where Pt is the transverse momentum of the backward proton and we used the light cone 
wave function normalized by 

dacfpT \4'{a,pT)\ = 1. (6) 



The neutron structure function F2n is (we shall suppose) known. The test we propose 
is to measure the cross section for backward nucleon production at a variety of x and a 
and examine the ratio 

Here, (Tmeas is the measured differential cross section and K is the factor 

A- = i^^=^(l + (l-.f). (8) 

The signature of a two nucleon correlation model is that this ratio is independent 
of X for any fixed a and px- Of course, how useful this signature is depends upon how 
different we may expect the result to be for a 6q cluster. This we shall see in the next 
section. 

2.3 Backward nucleons from 6q clusters 

For the case of electromagnetic scattering from the 6q state, we have basically the con- 
volution of F2 with the fragmentation functions of the five (or more, in general) quark 
residuum. Since the quarks in the residuum depend on which flavor quark was struck, 
we must write 

Here V^ is the distribution function for a valence quark in a six quark cluster, the sum 
is over quark flavors, and -Dp/5g is the fragmentation function for the 5q residuum, i.e., 
the probability density per unit z and Pt for finding a proton coming from the 5q cluster. 
It is tacit that the correct 5q cluster, either u^d^ or u^(P, is chosen. Argument z is the 
fraction of the residuum's light-cone longitudinal momentum that goes into the proton. 



the factor (2 — x)~^ comes because Dp/^q is probability per unit z in its definition and we 
quote the differential cross section per unit a. The formula is written for high backward 
proton momentum, where we can expect the 5q residuum and hence the actual 6q initial 
state to dominate. 

Neither Fg nor Dp/^g can be said to be known. However, since a large fraction of the 
short range part of the baryon number two state may be in a 6-quark cluster, we should 
make the best guess as to what these functions might be and see how large a difference it 
could make experimentally to have significant 6q cluster contributions, at least in given 
regions of phase space. 

Estimates of F2 in a model where the nuclei are treated as containing some fraction 
6q clusters have been given by Lassila and Sukhatme [§. They chose their quark distri- 
butions beginning with quark counting rules and then fine tuned with physical logic to 
describe the EMC data. For completeness, the three parameterizations they present are 
recorded in the Appendix. 

The fragmentation function is even less well known since there is no complete body 
of data to check it against. The counting rules suggest a cubic dependence, as (unnor- 
malized) 

Dp/gs{z) ex (1-;^)=^ (11) 

for z ^ 1 and for zero pt ■ We shall use this form, although we keep in mind the possibility 
that higher order contributions or renormalization group considerations could somewhat 
alter the power, as they do in many parameterizations of the quark distributions in 
nucleons. 

If the high momentum backward protons come from a 6-quark cluster, then ameas = 
aQq and the experimental ratio Ri should be given by 

o _ 0(6) _ F2\x) Dp/^q{z,pT) , . 

'^'-'^' - (2-x)(2-a)F,„(0- ^''^ 
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There is no reason for R\ to be independent of x for fixed a and px ■ We plot this -Ri 
in Figs, m and |^ for px = and specified a. Some old and simple quark distributions of 
Carlson and Havens [^ were used to get F2„ in Fig. [l| and the CTEQIL |ll| distributions 



were used in Fig. |^. The difference between what is seen and the horizontal line expected 
from a pure 2N correlation model is not negligible. That the ratio goes to a finite value 
as we reach the kinematic limit x = 2 — a in Fig. |l| has to do with the fact that both the 
5q fragmentation function and the dominant quarks in the nucleon approach their end 
points cubically in our models. A dive to zero or a flight to infinity is not precluded in 
real life, and the latter is seen in Fig. |^. 

3 Commentary 

3.1 Potential dominance of 6q cluster 

We expect that a correct description of the deuteron would have a neutron and a proton 
treated as in the 2N model when they are far apart. As they get closer, QCD processes 



such as gluon recombination |T2[ will surely occur and affect first the ocean parton distri- 
butions. It is something of a simplification to think of a deuteron as just a combination 
with a large fraction pure 2N state with a small fraction (perhaps 5%, from wave function 
overlap estimates [|T^) 6q state added in. 

However, we emphasize that while a 6q cluster may be a small part of the deuteron 
overall, it could be a large fraction of the short range part of the deuteron wave function. 
The deviation from what is expected in a pure 2N correlation could be as large as is 
shown in our Figures. The situation is not like the EMC effect where the effects of the 
6q cluster are diluted by the mostly ordinary collection of nucleons in the target. Here 
we can select events in a phase space region to enhance 6q cluster effects. The backward 
proton is a tag that emphasizes the 6q cluster and — if it is there at all — it will dominate 



the cross section for large enough backward proton momentum. 

It is necessary to have data at fairly high backward momenta. No one doubts that 
at low relative momentum or long distances the deuteron is a neutron plus a proton. 
How high is needed? We suggest 300 MeV/c is a good starting point, based on existing 
backward proton data from deep inelastic neutrino scattering and studies of the backward 
proton spectrum in that process using 6q cluster ideas. 

3.2 Falling {x) with increasing a 

Let us point out a piece of kinematics. From momentum conservation one has < x < 
(2 — a). Hence, unless the x distribution has a bizarre shape, one expects that {x)a — the 
average value of x at fixed a — decreases as a increases. This was pointed out in Ref. 
[^ in the context of the 2N model, and was initially suggested as a test of that model. 
However, the result should be produced by any model, so finding the trend in the data 
is not startling. 

The two-nucleon correlation model does give a specific result that {x)a falls to zero 
linearly as 2 — a. In contrast, the six-quark cluster model may or my not fall quite 
linearly. It depends on the specific implementation of the model. 

The rest of this subsection attempts to show why the two-nucleon result for {x)a is 
independent of the internal details of that model, but similar manipulations for other 
models do not lead to such definite results. It has to do with the way the cross section 
factors. 

For the two-nucleon model, the structure of the formula for the cross section differ- 
ential in X and a is. 



where ^ is defined earlier. An elementary calculation gives 



Jq "" dxx P{x,a) 
jQ-''dxP{x,a) 



{^U = Z2-:: :' ' : = (2 - «) (0 , (14) 



where {C,) is independent of a. We can easily turn this into 

7^ = (2 -a). (15) 

The sort of result that can be derived in the corresponding way for the six-quark 
cluster appears less useful. We envision one quark being struck and driven forward, and 
the residuum that remains "fragments" (or recombines) into a nucleon, that often goes 
backward, plus other stuff. The differential cross section is structurally 

P{x,a) = - — — = ^^^-^-D{z= ). (16) 

ax da 2 — X 2 — x 

The factor g{x) is for the quark knockout and D{z) is the fragmentation function. The 
first is a function of x since it involves distribution functions of quarks in the six-quark 
cluster, and the six-quark cluster is standing still in the lab. Argument z is defined 
earlier. The neatly derivable result is for average a at fixed x: 

, , in~^ da a Fix, a) , , , , , , 

(a). = V- , p/ '/ = (2 - x){z) , (17) 

Jq dar[x,a) 



or. 



(2-a;). (1^ 



{a)x=i 
This result for the six-quark cluster contribution requires averaging over all a, and 

while the optimist may expect the six-quark contributions to dominate at high a, no-one 
expects them to do so at moderate a. So, this result seems untestable. 

For (x)q, details would have to be worked out for each special case. However, remem- 
bering the general result that {x)a decreases with increasing a and is zero kinematically 
when a = 2 makes it likely that one will get something like the 2N result. 



3.3 PossiUUties at CEBAF 

Can CEBAF with a 4 or 6 GeV beam be useful? We believe yes for the 6 GeV beam. The 
question centers around how much backward momentum is possible with an incoming 
electron of this energy, and how much data we can get in the scaling region. 

3.3.1 Limits on backward proton momentum 

The maximum directly backward proton momentum when an electron scatters from a 
deuteron is (3/4)m7v, or 704 MeV, but this is for the case of infinite incoming energy. 
If the energy is finite, the magnitude of the maximum backward momentum is reduced. 
For example, with E = A GeV and Q^ = 1 GeV, the maximum directly backward proton 
momentum is 

p' < mOniN = 564: MeV, (19) 

which is still a decent backward momentum. (For a 6 GeV electron beam and the same 
Q^, the maximum p^ is 609 MeV.) 

These results were obtained with the help of 



aruN < - [md + iy- \/iy'^ + Q^] {l 



and also, for no transverse momentum. 



4m^ 



\ 2mdiy -Q"^ + m^ \ 



2 ' 



(20) 



p = TUn 



a^-l 



2a 



(21) 



For V — > cxD, we recover the limits cited, for example, in [Q]. Now at finite energy and 
for a given Q^ we maximize the limit on p^ by maximizing v. We have 



V = E- 



AEsm^{e/2) - AE' 



(22) 



For Q^ = 1 GeV^ and E = A GeV, we get a < 1.77 and p^ as quoted above. 
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Much of the hinitation actually comes because fixing Q^ for a given incoming energy 
puts a lower limit on Bjorken x. Since this is also the momentum fraction of the struck 
quark in the lab frame, it means that the struck quark is not moving forward as fast as 
possible, and the residuum is then not moving backwards as fast as possible. For the 
4 GeV beam and the present Q^ limit, Xmin = 0.14. If we went to infinite energy, but 
maintained this value of x, a would still be limited by a < 2 — Xmin = 1.86. 

3.3.2 The scaling window 

Using the formulas with the scaling result for F2„ requires that we be in the scaling region. 
This will squeeze the allowed range of x to be narrower than the kinematic limits. 

Scaling requires, at a minimum, that Q^ be above 1 GeV^ and that W (the photon 
plus single target nucleon cm. energy) be above 2 GeV, out of the resonance region. 

The lower limit on Q^ sets a lower limit to x, since for a given energy there is a 
maximum energy transfer u possible. The maximum u comes for backward electron 
scattering and leads to u < E — Q^/AE or 



AEQ 



2 






J mini- 



The lower limit on W sets an upper limit on x. Applying the limit to the final state 
that comes from striking the neutron in the 2N correlation model gives 

W^ = {Pn + q? = ml + 2mNv{2-a)-Q^ > W^,^, (24) 

for a given a of the backward proton and letting p^ = m%. This leads to a limit, also 
reached for 180° scattering, that 

_ 2E (2m^E(2-a)-(iy^,„-m^)) 

The curves giving Xmax and Xmin are shown in Fig. ^for a = 1.4. The scaling region 
is the region between the two curves. 
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For a = 1.4, the scaling region includes only a short span of x for electron energy 
E = 4 GeV, but the span increases greatly for E = 6 or 8 GeV. The span for 8 GeV is 
over half of the maximum possible at any energy, and not too much less than could be 
got at 15 GeV. For increasing a, which corresponds to increasing velocity of the struck 
neutron away from the photon, larger energy in the laboratory is needed to reach the 
scaling region. 

3.3.3 Outside the scaling region 

Do we need to be in the scaling region? Our formulas for the 2N correlation are simplest 
to evaluate there, and we don't know how to guess at forms for production off the 6q 
component, so our comparison case is gone. But, from the 2N viewpoint, the purpose of 
scattering off the neutron is to free the proton — nothing more. If all we want is a yes/no on 
the 2N correlation model, we could take a measured cross section for producing backward 
protons and divide by a cross section for scattering off a neutron at the correct values of 
the incoming variables, and see if the result depends only upon a and pt- 

The disadvantage of doing this may be more practical. Driving the struck particles 
forward forcefully separates them greatly in momentum space from the backward proton 
and reduces the final state interactions, which we have neglected in our discussion. As 
we leave the scaling region, it can mean that the energy transferred to the forward 
moving particles is low enough that we need to rethink our attitude towards final state 
interactions. 

3.4 Possibilities with muons at higher energy 

At Fermilab, for example, experiment E665 scattered 490 GeV muons from deuterium 
and xenon targets in a streamer chamber so as to be able to observe backward charges. 
A later version of experiment E665 has capability to observe backward neutrons, and 
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uses carbon, calcium, and lead as its heavier targets. The x dependence of the ratio 
Ri remains a good observable. There is a clear advantage in having higher energy as 
one does not have to think about the "scaling window." It is virtually the full possible 
kinematic range. 

Also there is some rate advantage in measuring backward neutrons instead of protons. 
For the 6q model, one can work out that if the 6q valence configuration dominates, as it 
should at high backward nucleon momenta, then the ratio of neutron rate to proton rate 
should be 3/2 and be independent of x. (One uses a combinatoric argument to get ratios 
of proton and neutron production from u^d? and u^d^ residua, and some discussion of 
this appears in |^.) The result for the 2N model is different, since one has 



0-2Ar(/id -^ finX) _ F2p(0 



(26) 



SO that the ratio should vary from about 1 at low ^ to about 4 at high C,- Thus back- 
ward neutrons are always produced more copiously, and this is an interesting additional 
observable if one has the capability of detecting both flavors of nucleons. 

We might also note that study of backward nucleons from a deuteron target is a study 
of the "target fragmentation region" and is best and most easily carried out in the rest 
frame of the target, i.e., the lab. Data presented in the photon-target cm. is not equally 
useful. 

4 Conclusion 

Study of the deuteron should be pursued since deuterons are our chief source of infor- 
mation about neutrons and we should understand this source. Further, the behavior of 
the deuteron state at short range gives information about the short range dynamics of 
strong interaction QCD. 

We have suggested a measurement to learn what the short range wave function of 
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the deuteron is. Namely, examine the shape of the measured differential cross section 
for electroproduction of backward protons or neutrons from a deuteron target, and take 
its ratio to what would be expected for deep inelastic scattering from a free neutron or 
proton, respectively. A high momentum backward nucleon acts as a tag isolating events 
where the initial material in the deuteron was tightly bunched. The x dependence or lack 
of X dependence of the ratio is a signal that is distinct for the extreme cases of a pure 2N 
or pure 6q cluster. We have presented simple model estimates of the size of effects that 
may be seen, showing that factors of two differences from maximum to minimum may 
be expected in the 6q case, whereas no maximum to minimum difference is expected in 
the 2N case. 
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A Appendix: Distributions for six-quark clusters 

We apply the notation |6g) or 6q to label the situation when the neutron and proton 
are melded and lose their individual identity. This notation emphasizes the fact that 
standard QCD quark parton model considerations should be applied to this interest- 
ing multiquark object. Therefore, for a generic \nq) state the sea, valence, and gluon 
distribution functions (times z) are written 



Gn{z) = Cr^il-zY'^, (27) 

where z is the fraction of the total cluster momentum. The coefficients and powers are 
determined in ||^ by appealing to standard normalization and momentum conservation 
considerations along with input information from experimental study of the n = 2 (pion) 
and n = 3 (nucleon) situations. As a result, 3 cases were developed to illustrate the 
sensitivity to small changes in the power of (1 — z), i.e., (a6,&6) = (11)9) for case A; 
(11, 10) for case B; and (13, 10) for case C 
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Figure 1: A putative ratio Ri assuming the measured cross section for backward protons 
in deep inelastic scattering is dominated by 6-quark configurations and using the LS 
model [9] for the distribution functions of the 6-quark cluster. In general, Ri is the 
measured cross section (sans some kinematic factors) divided by the neutron structure 
function F2„. If backward proton production were dominated by 2-nucleon correlations, 
the plot would be flat. The plot is for fixed a = 1.4 (momentum 322 MeV for directly 
backward protons) and uses CH [10] distribution functions to obtain F2„. The heavy curve 
is from LS parametrization A, the normal curve is from parametrization B, and the dotted 
curve is from parametrization C The vertical units are arbitrary as the fragmentation 
function -Dp/sg is not normalized. 
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Figure 2: Like Fig. 1 except that F2n is gotten from the CTEQ [11] distribution func- 
tions, specifically from the set CTEQIL for Q^ = A GeV^. The heavy curve is from LS 
parametrization A, the normal curve is from parametrization B, and the dotted curve 
is from parametrization C The curves turn up as x ^ 0.6 {^ -^ 1) because the CTEQ 
distribution functions all approach zero at the upper limit faster than (1 — ^)^, and the 
fragmentation function we use goes to zero as (1 — 2;)^. 



Figure 3: The scaling window for a = 1.4. Values of x between the two curves can be 
reached in the scaling region for a given incoming electron energy E, given in GeV above. 
The lower curve is set by the requirement that Q^ > 1 GeV^ and the upper curve by the 
requirement that W >2 GeV. The curves begin at E = 3.72 GeV for this a. 
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